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ABSTRACT 


A fast Fourier transform prograra for the generation 
and reconstruction of Fourier Transform Holograias .is 
developed Based on the algorithm suggested By Cooley & ■■ ■■■ •• 
Tukey .The theory underlying the Hologrsiii generation .& 
reconstruction and the algorithm for the fant Fourier 
transform are discussed. The savings in time and storage over 
the direct computational method are distinctly brought out. 
These savings are quite significant for large arrays. Certain 

applications of the FFT program in antenna pattern synthesis 
convolution integration etc. are also discussed. 
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1. INTRODUCTION 

Hologram techniques constitute a unique method of 
storing the three dimensional information (cunplitude as well 
as phase) so that it can be recall edat will as a real image 
that can be inspected in detail. Hologram photography may 
help in solving complex problems in image formation such as 
character recognition, aberration balancing of lenses. 
Practical embodiments of these hologram applications have 
yet to be accomplished, since these are based on the use of 
coherent illmmination.The feasibility studies can be carried 
out by simulating the optical process on a digital computer. 
The various steps involved in the generation & . ^ 

reconstruction of the Fourier Transform Hologram on a . 

digital computer are (a) generation of redundant data points 
by interpolation (b) computation of the Fourier transform of 
the modified data (c) addition of a uniform light pattern 
with linearly increasing phase to the trsnsforn (d) computing 
the intensity of the complex signal to get the Fourier 
Transform Hologram (e) inverse Fourier transforming the 
hologram to get the actual data without the redundancy 
introduced. 

Since the amount of data to be proeesced in image 
formation is usually prohibitively large from the viewpoint 
of computer storage, this calls for an efficient use of 
computer time & memory. The fast Fourier transform algorithm 
suggested by Cooley & Txikey^is one such method of efficient 
and economic computer usage. A Fortran I¥ program based on 
this algorithm has been written and debugged. The theory 
underlying Fourier Transform Holography, hoi ograi'ii generation 
and reconstruction technique, the details of the algorithm, 
the program amd its application^in antenna pattern synthesis 
convolution integration etc. are presented here. 



2, GEJ^EfiATIOK AiMD' JiECOKSTRUCTION OF HOLOGRAMS 

Fhe uagnitude and phast: of a scattered wavefront can 
be recorded photographically by superimposing a coherent 
reference beam or background wave on the field striking the 
photographic plate. The recorded diffraction pattern, called a 
HOLOGRAM, bears little resemblence to the object, but contains 
most of the information required to reconstruct the object. 
The process has three essential paxts: 

(a) The defocussing or spatial frequency dispersion 
of the image. 

(b) The holograia recordii% process. 

(c) Reconstruction. 

2.i Dispersion Process ; 

Let a transparency h£,ving amplitude transmittance 
s=s^+s^ be placed at plane P^.Amono chromatin plane wave 
illuminates the plane and is modulated by the 
transparency, resulting in the light amplitude distriWtion 



Construction of Hologram 



At plcUiu of focus" image or Fresnel 

diffraction pattern of the transparency is formed. The image 
using the? Prosnel-Kirchoff diffraction integral becomes: 

X(x,y)=|^ f}[s^+s (v.,v)] e-^T dudr 

•nao ^ “* 


Making the small angle approximation. 


where, a 


where 


[22+(x-u)2+(y.v)2J^ 4 +^2?^ 

' ' ' 

X(x,y)=^_^|s^+s^(u,v)j e~^T +iy-^) ] audv 

constant phase factor, e”^ ^ has been dropped. 
X(x,y)=s^+s^(x,y)*f(x,y) 

f(x,y)=^ e-^f 


X(x,y)=s^+ j|^S^(u,v) .F(u,v| . 

where 

S^(u,v)=jfs^(x,y)] 2 ;^ 2 2, 

and F(u,v)=5£f (x,y^ = e^ .a ^ 

2.2 Recording: 

Since the photographic plate records only the 
absolute magnitude, therefore, the recorded hologram 
represents the function: 

XX*=s^^+s^(s^*|:)+s^(Sj,*f)*+ js^^fj^ (-1, 

2.3 Reconstruction : 

The reconstruction is carried out by illtminating 
the hologram with coherent light. Considering only the real 
image term of Eq. ( i) , s^(s^*f)*, when this term is operated 
on by the dispersive filter, this becomes; 

Sb(s^*f)**f=3^[s^.Sr*(-u,-v)F*(-u,-v) .F(u, v)] 

: = -^b-^r' . 



This teriii along with the bias term.s, ^.coiatitutes 

D 

the reconstructed itiSge. 



The object at plane is iiaaged at plane P^ and the 
hologram is recorded at some plane P2 between P^ & P^.When 
the hologrciu is produced, the original object is removed and 
the hologram is reinserted at P2.The waves generated at P2 
are replicas of those that had impinged at P2 when the 
holograifl was made. These waves continue to P„ and form an 
image . 

If the phase function has the property, 

F(u,v)=P (-u,-v),the light distribution at the hologram 
plane will be real and in construction there will be no 
interference from the twin- image. 

The other method to eliminate twin-image 



interference is to put the dispersed signal on a carrier as 

shown in the Fig, above. The dispersed signal can be 

represented by: j| x 

X=s, +s^*f+s_o ^ 

u r c 



where 


I =_£ 

'c 

the recorded signal can be represented by: 


XX 


"^^b ■^2s^.s^.cos^^x+ |s^*f j^+s^[s^*f+(s^*f)*] 
+Sc[(Sr*f)e ° +(s^*f) e J 



If s^*f has bandwidth, W, the requirement that the 
spectra are non-overlapping is that 

2.4 Fourier Transform Holography ; 

For attaining high resolutions in wavefront 
reconstruction imaging, Fourier Transforia Holography is 
employed. The recording procedure for such holography is 
shown below: » u.u 


f hctjq» 



CoWiwvuter AoSC>t-bjyn 



Let the reference signal be: 
where 

0=bu+cv 

b=offset along the x-axis 
c=offset along the y-axis 

The complex amplitude at the focal plane of (by Pourier 
transforairtion) is given by: 

A(u,v)=H(u,v)+P(u,v) 

where 

H(u,v)= 1 h 1 J[h(x,y)] 

P(u,v)= JAqI 

The intensity recorded oii the. plate is given by: 

I(u, v)=A(u, v) .A*(u,v) 

= \E\^+ [Aq|^+ iHl-lAJ 

The reconstruction is done by Fourier trcinsformation of the 
recorded hologram, which gives: 

R(x,y)=reconstructed image 

= |H| )Aq| i-+ |A^.h.^(x+b,y-)-c) 

+ |Aj.h.5'(x-b,y-c) 

For negligible interference from the twin- image, 

b,c^|¥ 

where WxW is the width of the original object. 



3. COMPUTATIONAL STEPS INVOLVED IN THE GENERATION AND 

RECONSTRUCTION OF FOURIER TRANSFOM HOLOGRANIS 

To gen.fc;ra.tc.' & reconstruct the Fourier— Trausforra 
Hologram on a digital computer, the following steps are 
required: 

(a) If the data has NixN2 samples, redundancy is 
introduced by interpolation to generate 4NixN2 
samples since the hologram requires a minimum 
of four times as much band width as the signal 
to avoid twin-ima.ge interference. 

(b) This delta is Fourier transformed to get four 
times the number of Fourier coefficients 
required for the faithful reproduction of the 
original signal. 

(c) The carrier signal having a constant amplitude 
and linearly increasing phase is added to the 
trai'isformed array and the amplitude of the 
complex signal is computed to get the hologram 

(d) To reconstruct the original signal, the 
hologram thus genereited is Fourier inverse 
transformed. The inverse transform thus 
obcained contains the original signal, the 
twin- image, the auto-correlation function of 
the signal. Since the auto-correlation function 
of the signal requires twice as much band 
width as the signal.' and the twin- image 
requires as much band width as the signal ,the 
total (minimum) band width needed for non 
overlapping spectra of the signal, the twin 
image, and the auto-correlation function is 
four times as large as the signal bend width. 
Thus the redundant samples introduced mininize 
interference, and hence, the reproduced signal 



contrans the actual time signal without any . 
redundancy what so ever. 

Since the rmount of data to be handled is 
prohibitively large from the point of view of computer 
storage, this calles for an efficient use of computer time 
and storage. A multi-access, minimum time scheme is ideally 
suited for hcindling such large data arrays. The fast Fourier 
transform algorithm suggested by Cooley & Tukey is one such 
scheme. The theory underlying, this algorithm is explained in 
the following sections. 



4. FOUHIEfi THAWSFOM REPRESENTATIOI OF DISCRETIZED 

PERIODIC SIGNALS 

liot A(t) be a band limited sampled function of 
periodicity Tj-Sciiiiplwd at intervals of At, the Fourier 
transform and its inverse are then given by: 

X(J)= ^A(K)¥'^^ for J=o,i,2, /..^N-i 


where 


A(K)= 2 

T=o 

t=K.At 


for K=o,i,2,... ,N-i 


T=N.At 

¥=exp(27vj/N) 

A derivation of this expression is given in 
Append ix-B 

In the two dimensional case, the transform and its 
inverse can be written as; 


Na.-i w-i 


X(J.,Jp)= ZZ A(K.,K2)¥i 

X C. Iv't * X C- 


• Ji.Ki \,^J2.K2 


^ Ulso klsO 

for Ji=o,i,2, ,Ni— i 

J2=o, i,2, ...... . ,N2-i 

Mi-1 Nt-J 

-Ji.Ki „o-J2.K2 

A(K.,Kp)= 2J Z *^2 


for Ki=o, i, 2, ....... ,Ni-i 

K2=o,i,2,... ,N2-i 

where 

¥i=exp(2}lo/Ki) 

¥2=exp(2J'io/N2) 



5. C0MPUTATI01\iAL TIME INVOLVED IN DIRECT CALCULATION OF 

FOURIER TRANSFORM 

Thti number of operations in the one dimensional case 

2 

IS N , since from the very definition, N complex 
multiplications as well as additions are to be carried out 
for each point. 

Simila.rly, in the two dimensional case there are 
Ni.N2 mult iplicacions as well as additions for aech point 
and there a, re Ni.N2 points. Therefore, total number of 
operations to be carried out is given by: 

(Ni.N2)^=N^ 

where 

N=Ni.N2 is the total number of data points. 

Thus, in general, the total momber of operations in 
the direct calculation of Fourier transform equals the 
square of the total number of data points under 
consideration. 

12 

For N=2 ,for example, the fatal number of 

24 7 

computations needed equals 2 1.7xio .Since each 

calculation requires about 20 sec. on IBM 7044, the 

computation time required is about 340 secs. 

Presented in the followirng section is a. method based 

on Cooley & Tukey algorithm which requires Nlog 2 N 

computations, thereby, reducing the computation time 

iP ^ 

considerably. For example, for N=2 ,the time required is on 
only one sec. as against 340 secs, required by direct 
computation 



6. aOOIEY TUKEY ALGORITHM 
6.i One DiiaensioPLal Csise ; 

The Fourier transform in the one dimensional case 
can be represented as: 


N-i 

X(J)= 2'a(K)W‘^^ for J=o, i, 2, . N-i 

where 

A’ are the given complex numbers, & W=exp(2Aj/N) 

M ' ■ ' • 

If K=2 =Number of data points 

J can be written in the binary form as: 



where ‘^o’'^i’ ’’^'^M-i either 0 or I. 

Writing in the same way, 

K=%.j,2“-U 

The Fourier transform can be represented in. the ■ 
following form: 

1 1 

X( J^-i^ * • • * • • Z "'Z • ’Kq) • , 

* ^M-i O 

Since, =exp( 2xj )=i 

TV pM— i / j o^”^4. . j-T W qM— i T V’ 

■ W • ■ ■ =W . . . • =w 

The innermost sum in (S) depends entirely on J^, 

Letting Aj_(JQ,Kjy._2, . .. ..,Kq) denote this sum, 

and defining A 2 (JQ,Jj_,%_ 3 »- — -,‘^K^) as the neict' iniiermost 

sum over Kj ^_2 6n,we obtain the algorithm: 

■ 4 • J K - 2^~^ 

Ai(Jo»Kj^_2» • * ‘ ’^o^" ^^%-i^‘ * * ' * ’^0^'^ 

. . 1 

^2^'^o’'^i^^M-5' * * * ’^o^~ ^ ■^i^'^o’%-2^ ‘ ’ ‘ 


c&) 





JL 

*" 

W ^ ^ ° 

since, n, ■ 

’ pM-i 

W =-i 

~‘^l-i^'^o’ • ’'^L-2^-^^M-L-i’ “ “ ' * 

Jj-i 

( — i ) '^L— ^ '^L— 2^"^^ ^ * "^ * * * * 

(J, 22^-2.... ^+Jj2«-i 
^ L-2 o 

for J^ ,-=o,i and all values of J^,..., 

Jj*“ 1 0 

JL_2.an«i ^M-L-i’ »^o* 

Calculating for L=i,2, ..,M yields 

X(J]^_i,.....,Jo)=AM(Jo^ »‘^M-i^ 

The implication of (g) is that the bit configuration 
of the index of Ajyj is to be put in the reversed order in 
order to obtain the index of the corresponding X 

At successive stages the exponent of W represented 
in binary are: 

M— i_ 

K]vj_i(Jo>o, ,o)=(J^2 + 0 +. . . . .+o) .Kjyj_^ 

^M-2(J.,J^,o,..,o)=(J.2^"Vj^2^”^)Kj^^2- 


CV) 


cB) 


M-i 




,+J 


Fi-L 


0 


)K. 


M-L 


TK* 

Thus W has been factored as: 
^JK_^JK mod N 



K (j , 0 ,. 
M-4i. 0 


K (J 
W 


»,0) K (J ,J ,0...,0) 
M-2 1 Q. 

.W 

. ,J ,0,.0) K (J 
0 0 M-1 

...w 


J ) 

0 


Tlie t reer-giaph. represerctatiO'M aorrespands Ug; the, 
JK 

faetoririg of W as explained above. 

Immple 

lert- us coiosider' a one dimensional exfimple using 
8 daia points ^ 

ai'ncsp A ( J ,K - ^ A(K - , . . . . ,2; ) . 

I 0 M-2 0 M-l M-i 0 

[J .K .2 i 
0 M-1 

W 


Pon 1T=8 

M-3 J 

0 

A (J ,K ,K )=A(0,K: ,K ,K ) 

i. 0 1 0 ""10 1 0 

for J =0,1 and all o-tiier values of L ,2 , 

0 1 0 
Subst ituting all pofi&ible values of J ,1 ,K ,we obtains 

0 10 

JL;4000)=A(000 )+A(100) 

1 

A (I00)^(000)-il(l00) 

1 


A (010)=ii(010)4A(llO) 
1 

A (1!10)=A(010)-A(110) 
1 


A (001)=A(001)+A{101) 
1 

A (101)=A(001)-A(101) 
1 


C9J 



A (01l)iA(0ll)+A(lll) 

1 - 

♦ I . 

2nd Summation s 4 (J ». 2 +j )£ .2 

^ 1 0 1 

A (J" , J ,K )= ^ A (J ,e; ,k ).w 
_ 2 0 1 0 1 0 1 0 
Olc’ul'^'fring for J =0,1 & all other values of J >K , we obtains 

1 0 0 
A (000}=A (OOO)-fA (010) 

2 11 

A (010)=A (OOO)-ii (OlO) 

2 11 


A (00l)=^i. (OOl)fii (Oil) 

2 1 1 

(011)=A (001 )-A (Oil) 

2 1 1 

2 

A (100)=A (100 )iA (110). W 
2 11 2 
A (110)=;A (lOO)-A (Il0).¥‘ 
2 11 


2 , 

A (101)=A (101)4A (lll).W 
2 1 1 • 2 

A (lll)ssA (103 )-A (lll).t 
2 3 '1 

• Once A ’s are computed', A ’s a,re no more needed. 

1 1-1 

The same ,s peace, which was previously occupied by A 's,can! 

•: L -1 

now be alioted to A ’ s .More over, adnce each appears only 

1 L-1 

twice in the calculation of the corresponding xi ’s, these two 

L 

A 's can be computed first and can now occupy the place of 
L 

A 's which have been exhausted thus reqxiiring only two- 
L-1 

more locations f or the complete storage of Al(K). 

5rd Summation i (J ^ +.J 2,f^J }h 

■ 1 p 1 , n n 



Substituting tbe v'^lues of A (k),we obtinins 

0.1 (O ) ( 0 )^A (1 ) B ) -Ml (3 ( 4 ) 4 ai (. 5 ) 4 -A ( 6 ( ’^ ) 

i.. (1 j=4.(0}-A(l)4A(2)-A(3)-fA(4)-A(5)Vi(6)-A...7) 

3 2 2 2 

(2)=;.k(0 -a( 2)-A(3)W -fA(4)tA(5)W 

3 2 

-*■ 1 ( 6 )-/ (vOw 

2 2 2 

A (3)=A(0)-ii(l)W -A(2)+oi(3)W +A(4)-A(5)W 
3 2 

^ l {6 )W 

2 3 

.. (4)=A(0)+A(l)W+.i(2)¥ +A(3)W -*.(4)^(5)W ^ 

^ -.i( 6 )¥ -A(7)W 

2 3 

A (.*i)-A(0)-A(l)W.^i(2)w -A(3)5ir -A(4)tA(5)W ^ 

^ ^H'6)W +Ji{7)W 

3 2 ^ 

A (6)-^(0)+A(l)i? -A(2)y -A(3)T?-A(4)-A(5)W 

^ t--i(6)¥ -A(7)W 

3 2 ^ 

A (7)=A(o)-A( 1 JW -A(2)W -A(3 )W-.ii.(4)+-i-i-(5)¥ 

3 2 

1-a(6)¥ +A(7)iv 

In tbe matrix notation, the iiario.us summations can 




x(ooo) 

x(i:.o) 

X(olo) 

X(ilO) 

x(ooi) 

X(lol) 

X(oll) 

x(lll). 


A (000) 

2 

1 

0 

A (001) 

2 

0 

1 

A (010) 

2 

1 

0 

A (Oil) 

2 

0 

1 

(100) 

2 

0 

0 

.i (101) 

2 

0 

0 

A (110) 

2 

0 

0 

A.^(lll) ^ 

0 

0 

T* ^ 

f- 

0 

A (000) 

3 

1 

w 

4 

A (001) 

3 

1 


A (OlO) 

3 

0 

0 

A (Oil) 

_ 3 

0 

0 

A (iOO) 
3 

0 

0 

A (1011 
3 

0 

0 

A (xlO) 
3 

0 

0 

[a (111) 

. L? 

0 


3 


0 

t 0 0 0 

0 

0 t 0 0 

4 

?? 0 0 0 

4 

0 W 0 0 

0 0 1 0 

0 0 0 1 

0 0 10 

0 0 0 1 


0 0 0 0 

0 0 0 0 

2 

1 0 0 
6 

1 W 0 0 

0 0 1 ^ 

5 

0 0 1 W 

0 0 0 0 

0 0 0 0 


0 

0 


Aiiooo) 

0 

0 


Aicooij 

0 

0 


AiColD) 

0 

0 


A IC0 11) 

2 

w 

0 


Attloo) 


2 



0 

w 


Ai(1^>i) 

6 

w 

0 



0 

6 

_ 


AixMl) 


0 0 


Ai.(0O]D) 

O' 0 


A (001) 

2 

0 0 


A (Olo) 

2 

0 0 


A (on) 

2 

0 0 


A (100) 

2 

0 0 

3 


A (101) 

2 

1 w 

7 


A (no) 

2 

1 


A (111) 



2 



The ¥ector J'x(a)J is merely with the hinary ■ 

arguaent flipped, therefore, 'by hit reversing the arguieent of 

vector I^Aj(n) one knows Oj^ectly where each component of the 

vector 1^. :)C(n) actually "bei^-ongs. The general fora of the 

factored matrices has heen developed by E.O. Brigham & 

2 

R.E. Morrow .This is discu=!Bed in Appendix-E. 

6.2 Two Dimensional Case ; 

Letting 

» M 

cl 

N^=2 be the number of points in the a.=i,2 
directions. Defining: 

Mi'=Mi 


M2'=Mi+M2 

we have -the ' indices Ji,J2 in the bit positions 0 to Mi’-i, 
Mi* to M2*-i, respectively, of an index J.Thus the bit 


configuration of J is described by the bit positions: 

M2' ,M2'-i, . . . . - ,Mi' ,]yii'-i, . . . . . ,0. 


,-r 2^^’ + T 




The bit configuration for K may similarly be written as: 


„„ .M2’-i „ A.V 

The two dimensional Fourier transform may thus be 

represented as follows: 

1^1 i 

X(J2,Ji)= A(K2,KP - 

wi ph2*^i. .J ^i’lKp 


wher« 


W =exp ( 2 Aj /M ^ ) , a=I , 2 

‘3i ■ ■ " Ci» ■ 



The L th summation over KI may be written from (7) in 
the following form: 


/ c- x’i -4- ~ ju 5 » • y j « 



^I^^M2’-i’ * • * 

’’^Mi’ ’'^o’ * • 


i 



= ^ \-i^^M2 

' _i> • * »^]y[if 


‘tvi-C*» 

Wi ^ 

letting 

3«=p.2^^+(l.'2^^' 


where 

p-0,I,2,.... . 

. , (2^^“^Ti) 


q=0,I,2, 



r=0,I,2, 

.,(2“'^~^-i) 

Putting 

in the values for p,q,r,J,’ 



^)=Vi(j)+. 

where 

Q=J^2^“^+ 

1 

•'^‘^L-2 



...+Jo 

Solving 

for , i , we 

obtain 


(jn 


,r?,oMi-L> 


(^J 2^ 1) 2^1-1 

- r ♦ Jj” 1 

wTi 


\ ^ ^ ^ =^-i ^ j >+-^L-i ^ ) Wi 


02) 


Q.2^^”^ 

d5) 

-r ^Mi-L 

0^ 


since 



2ad SuiiLtatioa : 

Letting A’ denote the result of the first suLiraation 

’ =-/']y[i(K]Jj2_i» • • . rKni, Jo» . . . , Jjyji^i) 

M 2 -I 

X(J2,Ji)= ^ A'.W2'^^^^ 
la-o 

11 

W 2 J 2 K 2 

Perfo ruing the siiuuation over Kjyj2'-L denoting it hy A£, 
we obtain: 

^'L ^ “^Mi ' • * ' ’ '^Mi+L-i ’ ^M2 ' -L-i ’ * * ’ ^Mi ’ *^0 ’ * * » '^Mi-i ^ 


05 ) 


Letting 


Where 


^L-i^'^Mi* " ’‘^Mi+L-2’^M2_^’ ‘ ' ’^Mi’'^o’ * * »'^Mi-i^ 


00 


W2 


/j pL-i j w pM 2 ->L 

'‘‘^Mi+L-i^ ■^•* ‘^Mi''^M 2 ’-L 


7 =p. 2 ^ 2 -^^^ 2 "l 2 --L-i^^ 


p=o 




Mi+L-2 


** T r)X)““2 J qIj’^3 , , T 

^"'^Mi+L- 2 ^ .■^'^Mi+L- 3 ^ ^ 

^M2-L+i 


Mi 


^"^^M 2 '-lH-i^ 


^ „Mi , , ^Mi-i , , T ■'- 

+ .. .. .+Kjy^)2 +(Jq 2 +*'+JMi-i 


Substituting for obtain the algorithm: 

MO» T ri o^ 2 -L 

A£(J)=A£_j_(j)+^£^i(j+ 2 ^^ "^)¥ 2 '^- 

— M 2 ^L 

/£( 5 + 2 ^ 2 ’'^)=A’_i(j)-A£_i(J+#^^ 


ov 

m 



where. 


corresponding to bit pGsitions:M2'-i', . 

9 = 0 . 1. 2 


. ,Ma’ 


corresponding to the bit positions: 

Mc’-i,,. .,.,MF,.'-L+i 

r=0,I,2, ,(2’^^-'"^-i), 

corresponding to the bit positions: 


Ma’-L-i, .. . . . ,1,0 


0 = 1 , 2 , 

corresponding to the two dinensions respectively 
6.3 Couputational Tine Involved In Cooley & Tukoy Algorithm : 

■fine involved in Cooley & Tukey algorithiii is best 
explained by reffering to thvj generalized factored matrices 
representation of the algorithm. As there are only two non 
zero entries in each row, and one of those two entries is 
unity, therefore, in each mtrix nultiplication the nuuber of 
couplex raultiplications and additions needed are 
respectively, N/2 & N, since in half of the terns, the conpl ex 


Laulti plication needed is already carried out in the other 
half of tile terns .Thus the total number 

of couplox multiplications and additions needed are, 

IVl 

respectively, ^N.M & K.M (N=2 ).The computation tine needed 
in Cooley & Tukey algorithn is thus roughly proportional to 
Nlog 2 N, where W is the total nunber of data points. 

Por N=2^^,the tine required on IBM 7044 works out to be 
nearly one second. 

Thus through the use of this algoritiin the tine 
roquireiient has been considerably reduced for large arra.ys 
(one Second as against 340 seconds for direct conputat ion 
on IBM 7044 for 4096 sejoples). 



7. PORTjS/iN IV PROGiiftM FOR IBM 7044 

7 . 1 Options : 

The program for Cooley & Tukey algorithm ha,s been 
written in Fortran language using the least sophistication 
level. The program can he made much more efficient by the 
use of MAP subroutines which offer the following distinct 
advantago'S over Fortran IV; 

(a) Compilation time is saved since Fortran IV is 
first converted into MAP and then processed. 

(h) A lot of time can be saved by the use of logical 
operations such as bit inversion, laultlplicotion 
uy <i , binary addition. These operf ,1ions are 
allowed only in MAP. 

(c) Computer storage can be nore eifectively utilized 
The Fortran program allows for calculating the 
Fourier transform as well as the inverse transform of one 
dinensional as well as two dimensional data(see nppendix-F 
fur Control cards). Print statements may be added to print 
out INTERMEDIA.TE location indices to identify the flow of 
the program at each summation stage. 

7.2 Capacity : 

Since IBM 7044 has 32K memory and since the program 
handles data points which are a power of 2 (2 , 2 , 

allowing for arrays AR(I) , AI(I) ,PR( ISC) ,FI( ISC) a total 
storage of I6K, which leaves I6K for program instructions, 
the Liaximum array size that can bo handled is thus 4096 
data points. The computer may more efficiently be used by 
one of the following techniques: 

(a) Simultjai-ous Fourier analysis of two sets of real data ; 

N-i 

X(J)=X’(J)+ 3 X''( J)= Z for J=0,I, . .,M 

|«:-o 

X’(J) & X'’(J) are real. 


where 



The coL.plex aKplitucles of 


X>(J)= Z A‘(K)¥'^^ 

N-d 

X"(J)= Z A"(K)¥'^^ 
are given by ( Appendix-C) 

A'(K)=4 A(K)+A*(K-K)] 

//'(K)=i A(K)-/*(N-K)] 

Thtsri;;fore, given two sets of real data, one can let 
these be the real and iiuaginary parts of the coiiiplex data 
and then the inversion formula can be used to obtain 
complex A' .These, used in a two dimensional of ('#^),will give 
the desired aiaplitude of {%) & (il) . 

( b) To Double the nuuber of data point s ; 

A Fourier analysis of a large number of data points 
12 

(larger than 2 ) can be performed by the use of external 

storage a.nd multiple references to the program. The procedure 
fur doubling the number of data points is as follows; 

Let N be the number of points and N/2 be the 

capacity of the core memory available. Let 

W=exp(2Ao/M) so that W^=exp(2Aj/-|-K) 

Then, in two separate passes with the program, one can 

ca.lculate the Fourier coefficients of the two N/2 Fourier 

transforms: u/ 2 _i 

X(2J)= Za'(K)(W^)'^^ 

'iZ for fcO,I,2,....iK-I 

X{2J+I) 2 A’'(K)(¥^)'^’‘ 
u=0 

The coefficients of_^ 

X(J)= ^ for J=iO,I,2, . . . ,N-I 

are given by (Appendix-D) 

A(K)=ifA‘(K)+A''(K)W^j 

for K=0, I,2, . . . ,4-1-1 
A(iN+K)=4-[A'(K)-A"(K)W”^j 



Repetition of the above procedure pernits successive 
doubling of the capacity of the program. 

7.3 Truncation Error ; 

To test the accuracy of the progran, several input 
wavefoms were transformed; The inverse transformation was 
perforued on the result and the r;n.s difference between the 
final and initial data was computed. The results obtained 


wer«^ as follows; 



N 

Type of date 

R.M.S Error 

64 

square wave 

M 

O 

1 

^ oo 

4096 

sinx/x 

lO" 



8. RESULTS OE HOLOGM GEHERaTIOU 

Two dimensional Eourier Transform Hologram 
wss generated' corresponding to a completely wMte 
transparency (encH sample Has unity amplitude) .Tiie Eourier 
ttrans-form of sucH a picture corresponds to an impulse 
at is clear' from the Eourier transform results. The 

reference signal consisted of a plane wave a minimm ten 
times brighter than the picture''s brightest, sample in 
order to keep the reconstruction error with in .This 
error is inversely proportional to the square of the 
brightness of the reference beam relative to the picture's 
bri^test signail samlpe, The relative figure of ten was 
aasrmied in tids particular exa.mple.The hologrR,m thus 

generated wfis inverse Eourier transformed to get the 

—6 

original aignal.An error of the order of 10 was obtained 
f or 8 p'Oints. 



D^TA POR TWO OIMPRSIONAL HOLOGRAM SYNTf^ESIS 


1 .00000000 

l.OoOOOOOO 

m 

1 .00000000 

1.00000000 



1 .00000000 

1.00000000 

l » O 00 000000 

1.00000000 



1.00000000 

1.00000000 

1.00000000 

1.00000000 



1 .00000000 

1.00000000 

1.00000000 

1.00000000 



l.OuOOOOOO 

1.00000000 

1.00000000 

1.00000000 



1,00000000 

1.00000000 

1.00000000 

l.OuOOOOOO 



1 .00000000 

1.00000000 

1.00000000 

1 . OoOOOOOO 



1.00000000 

I.OOOUOOOO 

1.00000000 

l.OOuOOOOO 




iOURIER TRANSIC.IM 
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0. 0000000 1 
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0.00000000 

0.0000000 

o.oooooooj 

0.000^00 00 

0,00000000 
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0.00000011 



9* GTHEH APPLIC/.TIOMS OF FFT PROGRAM PACKAGE 

The FFT prograia has also been employed to find the 
convolution integral and auto-correlation function^. Another 


important application of FFT program is in the synthesis of 
antenna patterns with equally spaced elements. This 
particular application of FFT is discussed below: 

For a general linear array with equally spaced 
elements, the relative amplitude of the radiated field 
intensity is given by: 

E=I+I.e^'^+ ,+I 

0 1 XL-1 

where 


and 


where 


are thw complex element currents, 
V=^*c’cos0+o( 

d=spacing between successive elements 
o^=progressive phase shift from left to right 
Eq.(2.2-) can be written as: 

n-^ 


E= 


Z=e 


dH' 


If we take discrete samples of E(J)' at regular 


intervals, then the phase associated with the J th sample is: 

~,J, J=0,I,2, ,n-I 

H 

The intensity for the Jth sample may be written in 
the following fora: 


where 


n-i 

E(J)= jT 


Eq.(A3) is the complex Fourier inverse of the 
function Ij^.The required comiplex current amplitudes can be 
found by the Fourier transform of the discrete function E(J) 



Tht3 CDLiplux current 8jj.pl itudos are thus given by; 

1^4 J E{j)vrJK 

3=0 

Thus, a. a;iven antenna pattern can be transforued into 
the (lonain, (iiscretized and Fourier transforued to get the 
complex current Giaplitudes,!^^, which produce the desired 
radia.tion pattern. The specific case of an end fire arra,y 
with 'X/2 spacing between successive elements is discussed 
below; 

For d= yz 

'f'=^COS0+oC 

For an end fire array: 

(X =- A 


1*= A(cos0-1) 

The simple pattern shown in Fig. 7-1 when transformed into the 
domain looks like Fig.3’2 

which is the Fourier series expansion of the function f(*v) . 

sin K ^/2 

This pattern has been synthesized for N=64 using — • , 

K=0,I,2, ,63, as the input data. 

Similarly, synthesis has been carried out for 64 


element two dimensional array using 


sin K'^/2 sin K”a/2 


K'A/2 


K"A/2’ 


K*=0fl,2, . .♦,7,K''=0,I,2, . ..■,7,as the two dimensional data. 
The results thus obtained have been 'plotted against the 
ideal patterns expected for 11=00. 





10, OOHOLUidlOdS 


One and two dimensinnal Poiirier Transform Hologram 

weire auccessfully generated & reconstruct ed on IM 7044 

ua±ng PKE program package. The results are presented in sec, 8. 

Synthesis technique for one- & two dimensional antenna patterns 

for equally spaced elements using PPT program package has 

also been successfully demonstrated . 

The value of the PPT program lies in the reduction 

of the computer time. An H-point transformation hy the direct 

2 

method requires a time nroportional to H where as the PPT 

program requires a time proportional to H. log N. The 

2 

approaomate time; saving factor is thus given bys 

H/log I 
2 

Other importent applications oi the FIT program 
package consists in finding the convolution integial,the 
auto-correlation function et-c.It has made -’-possible analysis 
in the frequency domain which was considerei impracticable' 
otheivrise.This may find enormous applications, in band width 
reduction schemes, character recognition, spatial filtering. 
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Programming Details 


APPEMDIX-A 


The summation over the first dimension is given 
by: 

~?EEEEI5E~ ^ > ^ i z- n z zzrz-r--" 

“ '^-i^%2-i»- •••’^Mi»'^o’ •»'^L-2»°’^Mi-I-i’*°*»^o^'^ 

^ '^L-i 

'^-i(^M2-i’ * * * ’^Mi*'^o’ * * ‘'^L-2V^’^Mi-L-i’ ' * *^o^ * 

(.Jj p.2^~2+.......+j 

Wi o 

The L th summation over the 2nd dimension is given 
by: ^ 

^*^Mi+l-2*'^MiH-L-i* %2~L-i* • " ^ *Sli *'^o*'^i* * * »'^Mi-i) 

> ' < => 

“ \-i^'^Mi’ ''^Mi+L-2*^’^M2-L-i" * * * »^Mi*‘^o*‘^iV‘ * * '^Mi-i^ 

+ (-i) \-i^'^Mi’**’'^Mi+L-i’^’%2-L-i^'**%i’'^o’V’ 

‘^Mi-i^* 

j-T 'l pM2— L 

^2^'^Mi+L-i*^ ■Mi''*^ 


The final summation over the two dimensions gives 
the Fourier transform with inverted bit configuration, i.e., 

^ ^^M2~i* • • • • • • • • *’^M i* ^F|i--i* • ‘ ‘ • ’ >^ 0 ^ 

=A^2^ *^Mi* • • * ^ • *'^M2-^i ^ ‘^o* ‘ • ’ * * 

From the last equation, it is clear that if the final 
transform has to be in the proper bit configuration, the 
summation at each stage has to be carried out with the^^ ^ 
required number of bits in the inverted configuration. This 
is shown by arrows in the above algorithm. The arrows point 
towards the bits requiring the lesst number of changes at 



each stage. This is illustrated below in the case of a 8x8 ’ 

picture plane requiring three bits in each direction, 



)+A(K^ ,K^ ,K^, 0 , Kq)"!' 

^ <1 ( 

) • A(K^ ,K^,K2 , i ,K^ , K^) 


).A^(K5,K^,K^,J^,0,Kq)+ 

<= ^ ( 

) • Ai(K^fK^,K^,J^,i,K^) 

(K^ f f f J ^ ^ 

).A2(K5,K^,K^,J^,J^,0)+ 

1 

).A2(Kj,K4!K^,jQ,J^,i) 


).A3(0/K4,K3,J^,J^,J2)+ 

i — - — ^ ^ ( 

).A^(i,K^,K^,jQ,J^,J2) 


) • Aj^ 

— ^ ‘ 

).Aj^(J^,i,K^,jQ,J^,J2)/ 

Aj(J^,J^,J5,Jq,J^, J2)=( 

) . A2( J Jq» Jj^> J2^'*’ 


The arrow on the top indicates the bit over which 
the sumciation is being performed. The problem is then 
reduced to calculating proper indices at each summation. 

let PP( I) denote the bit positions in the 2nd 
dimension over which the summation is yet, to be carried 
out, QQ( I) denote the bit positions in the 2nd dimension 
over which the summation has already been performed. 
Similarly,let R(I) denote the bit positions in the ist 
dimension over which the summation is yet to be carried out, 


Q(I) denote the bit positions in the ist dimension over 
which the summation has been carried out. The indices of the 
inverted bit configuration at any stage can then be 
represented as; 


and 


2 QQ(I).2' 
1-4 




,M2-I 


2 * 

I-* ^ ^ 


^ ^fi'(I) .2^'^ ■ 1 


whtrt L tli 
suiamatioa 


suumation in the first diinonsion and the Jth 
in the second dimension are being considered. 



APPEHDIX-B 

Ihe one dimensional Fourier transform pair is 
defined as: « 

A(t)4Xx(j«)ej“*d« 

If A{t>) is periodic with petiodicity,T, the transform 
pair becomes: ^ <x» 

. ,(t) 1 /X(3«)el-n(» -nu)^)d<J= ^ X(onw 


where 


•IK ^ 


r\~-» 


Letting 


where 


0)^=2 A/T 

If A(t) is sampled at intervals At 
t=KAt 

3n2,.J(At/.T 

A(KAt)= X(jnn>^)e ~f - 
N=T/At 

A(KAt)=2 

= ^ X(dnu> )¥^ 

"tvs U 


W=oj2*/N 

A{KAt) representsthe amplitude of the Kth sample 


and X(jnu»^) the complex amplitude of the nth Fourier 
coefficient, These will henceforth be abbreviated as: 
A(K) & X(n) ^ 

A(K)= Z' X(n)\/ 

For the periodic, sampled signals 
1 


,,nK 


X(n)4 


dt 




'T 'o 

M-l 


4 2A(K)e 


-jn2AKAt/T 


kiso 


ia-1 


'd[At 


2a(K)W 


riiK 





■ majc 


^ saiapling theorem, any function A(t) 

band limited to B cycles is completely described by the 

saa.ple values i/2B sec. apart. 
t=i/2B 

="/2\ax-^o =T/2n„ 

\ax=^/2At =iN 
Thus frequency band from -^N+i to is sufficient 
to completely describe the signal. 

N/Z 

A(K)= 2* X(n)W*.= Z x(n)W* 


/ nK 
X(n)W 

n*- N^+i 


Z X(n)¥^+ Z X(n)w’^ 

n? -N/a.41 

(m-N)K _ 

X X(m-N)W ■■ 

^ ,,nK N-l ^ 

^ ^ Za(K)W-<“-«K 

N:;! 

W •«.? . /..V . -nK 


2' X(n-K)w“^ 

nrN/a,4l 


J 2a(K)¥ 


^ NAt _ 
05>»fc.-vl U.=<1 


M-1 

Z 

n=»»/244 


.nK 


Substituting in the expression for A(K),we obtain 


Nil 


A(K)= 2' X(n)W*^ + Z X(n) 


,nK 


N-1 


nto 

N/2 


W 


,.,nK: 


n::i, 


Thus 


A(K)= 


2 ^ X(n)w“^ + ^ X(n)W 

M-1 „ 

Z X(n)W^ 

n*o 

N-i _ 

^ X(n)/^ & X(“^)=iAt 


M-l 




,-nK 


liso 

are Fourier transform pair for sampled periodic band limited 
functions. 



APPENDIX-C 


N'l 


X(J)= X'(J)+jX"(j)= 2 'a(e)w^k 

U.=o 

M*1 H-i 

= 2 ^A'(K)W'^^ + j 2 ^ A"(K)W‘^^ 

U-zo U .-0 

Whi M-d 

M'^)- Z = Tx’ (j)vr*^^+3 “T x"(j)w‘‘'^ 

- 3 ^ “^1 

M-A >4-i 

A(N-K)= ^ X(J)W“'^^^“^^= jKJjr"., 


-JK 




7^0 
H-l 

= Zx(J)l^' 

7=0 
rJ-t 


JK 


N-i 


= 2^' T X'’(J)W'^ 


because X’(J) & X''(J) are real. 

N'l N -1 

A*(M)= ^ X'(J)W"‘^^-3 2 X"(J)W"'^^ 
i|A(K)+A*(K-K)j= ^X'(J)W"'^^ =A'(K) 


N -1 


It O', 


i[A(K)-A*(K-K)| =3 2 X'’(J)¥"‘^^=3A’»(K) 

-* 3*^0 

A’(K)=i[A(K)+A*(N-K)] 

a»(k)=|j[a(k)-a*(n-k^ 



APPENDIX-D 

M-i M-1 

X(j)= Z A(K)wJK . 2 A(K)wJK, 2 A(K)wJ«' 

Letting fe *«/2 

m=K-iN 

N/i-l Wfc-1 

x(j)= 2 A(K)W'^^+ 2" 

U*** nv'j«=o 

M&-J N/a-i 

= ZA(K)r‘'+(-i)J TA(iH+lii)W'^“ 

) '"=« 

o T>r Wb-l 

X{2r)= 7 a(K)(w 2)^K+ 7A(iN+m)(w2)™ 

l/xA 

Nfe-l N/a-1 

= 7 [-^‘(K)+A(iN+K)] (W^)^^= 

Similarly, k=6 

X(2r+i^= W^[A(K)-A(iN+K)] 2 -^"(K){wV^ 

Compairing, we obtain? li=(> 

A'(K)=A(K)+A(i-N+K) 

A’'(K)= [A(K)-A(iM+K)] 
or j 

A(K)=i[A' (K)+A''(K)W“^] 

A^7W)=i[A’(K)-A''(K)V^“^j 



APPEND IX-E 

THeE-GHaPH eephesentation oe pactored matrices 

The generalized factored matrices can be obtained 
with the help of tree-graph representation. Taking asan- 
example, N=8, the factored matrices already developed in 
sec. 5 can be represented as a tree-graph as shown below: 



The sampled data A(K) are represented by the Q th 
vertical column of nodes on the left of the graph. The ist 
vertical column of nodes corresponds to vector Aj_(K),the 2nd 
column corresponds to A 2 (K) and the last column corresponds 
to A^(K)=X(K) .Each node is entered by a dashed and a solid 
line and within each node is an integer. The solid line 
brings a quantity from one of the previous nodes, multiplies 
the quantity by W^, where p is the integer in the circle, and 
the praduct is added to the quantity brought by the dashed 
line, For example, at node A^(ooo), 

Aj^(ooo)=A(ooo)+W°. A(ioo) 

A generalized method of establishing a graph of the 
form shown above for any N=2^ has been developed by Brigham 
and Morrow^. The generalized graph is formulated by the •. 'c 
following set of rules: 

X. Assume H sample values of the time function A(t) 

constitute the first vertical array of nodes,denot6d by A(K) 


where K=o,i,2, jM-i represent the sample times or 

address (locations) of of the function A(t) ;argument K is 
expressed as a binary number. This binary number determines 
the location of node A(K). 

The other arrays, Aj^(K) , are drawn successively to the 
right, and the nodes in each array are addressed as binary 
numbers as mentioned above. At each node, a number is written 
as explained by rule 2. 

2. The number in the circle of the Kth node in the L th 

array is found by; 

(a) writing the nuraber K in binary 

(b) scaling or sliding this number M-L places to the 

■ filling in the newly opened bit 
positions on the left by zeros. 

(c) reversing the order of bits 

For exauiple, refer ing to the graph for N=8,the node 4 in 
column 2 has the number determined below; 

Node 4 expressed in binary ; ioo 

Nmiber of right shifts required ; 3-2=i 

Binary number shifted one bit to 
the right & filled by zeros in ; oio 
the blanks created on the left 
Binary number after bit inversion ; oio=2 
Siuiilary , the 3rd node in the 3rd column has the following 
number in it; 

Node 3 expressed in binary ; oii 

NurAer of right shifts required: 3-5=o 

Binary liumber after shifting : oii ^ 

Binary number after bit inversion ; iio=6 ^ 



5* In array, node K (in binary form Kj^ . ,K^) 

has a solid line drawn to it from a node in the L-ith array. 
The address of tne node in the l-ith array is the same as 
node K, except that bit must be a one. The dashed line 

comes from a node ixi the L-i':'''.. array whose address is the 
aauie,but bit must be a zero. 

For exai.iplo,let 

M=5,L=2,and K=3=(oii) 

The solid line cones from ooi=3rd node of the ist • ■■ 

column. 

The dashed line comes from ooi=ist node of the ist 

column. 

4. The Fourier transform values are found in the last 

array, the location or argument of these transform values is 
flipped in binary order. For example, tne traxisform component 
located in ooi is actually the frequency component ioo, which 
is the bit reversal of the argument ooi. The address of the 


last srray must be bit reversed before the final transform 
values are obtained. 

The tree-graph merely corresponds to the factoring 

of W’^^.From Eq.( ),W'^^can be factored as: 

mod N 





¥ 





For M=3 

,,JK ,,JK mod 8 

W =2:W 



In the tree-graph the node A(K),the K^ time sample, 
is connected to any of the spectrum nodes A^(K)hyM lines. 



For tixaiaple, A( ooo) & A^(ooo) are connected only "by three 
dashed lines; these M lines correspond to factoring of 
The solid or dashed characteristic of the lines depend on 
the value of Kj^.If K^=I,the line is solid; for K^=0,it is' 
dashed. 

The Fourier transform values with the node 
end up in position Starting in the Q th array at 

position K 2 >K^,Kq, proceed by the solid and dashed lines to . 
location S®'^ contribution of node K 2 »Kj^,KQ 

to the transforued nude J 2 , J^-For example, start at locati 
location oii and head for ooi .to get the contribution of A(o 
A(oii) to the Fourier component A^(ioo). 

Array Position Power of W 

0 K2,K.,K^(oii) 

1 Jq,K. ,KQ( oii) J^,0,G(oe-)^0 

2 J.,J^,0(ooo)=0 

3 ,J2(ooi) J2,J.,JQ(ioo)=4 



iJutJLL GuOa,NAME K B OHRI ELECTRICAL 

SIBJOF 

f 

SIBFTO' MO DECK 

DIMENSION AR(6A) ,An64) ,FR(64) ,FI (6A) ,C( 64) ,AMP(6L) ,DIFF(64 

/'-iVv.';/-’: 

6 ^ , 0,( 6 UQQ { 6 )', R ( 6 ) , RR ( 64 ) , R I ( 64 ) 

P';' FORMAT (21 5) 

C .'-'.AKE M the LARGLR INTEGER 

700 IF(M-N2)701,702,702 

70X NTEf',R = Nl 

N1 = N2 
N2=NTEMf' 

H* 

PJ,I = 3. 14X59 

. XM = C.i+ALOG(XNl)/ALOGT2. ) 

XM2=0.1+ALOGCXN2) /AL06(2, ) 

'41 = XM1 
M2=XM2 

PRINT1,NX,N2 
FRINT1»MX,M2 

C DAT/. GLN. HAS STARTED. 

D02i=X»N 

ARr|)=l. 

A I ( I ) = 0 . : ■ 



Xlfl-'l 

RR{ I ) =10.*CCS(XI«PAI/i6. } 
Rid ) =10-.*SIN{X1*PAI/1'6. 3 

pRif.nioo 

. ^ , 



D08113 

D= ( AR { I ) ) «*2+ ( A I { I ) 


81i C(I)=SQRT(D) 

C DATA GE,J. !S OVER 

LCITECK=1 

uoo icc=o ; ■ 

V, y U being USE4 WIT'T N2 WHILE L IS BEING US#' W^ 

''J*0 


L = 1 



Z2=]. _ 

SS2 = 0 • * 

G TO TAKE CARE 06 ONE DIMENS IONAl' CASE 

3D9 IF{M2)S 08,308 ,307 

3 c 8 S ~ 0 • 

SS2=0. 

21 = !• 

^ 22? ' • 


307 lF(02-l)30U302»3Da 

302, DC303'!^1,J2 ' , 

•>03 QQ(I)~0» 



<=I-1 

ZZ=ZZ*QQ(I) 
S2«S2-H(2.^f*K)^^QQ( I ) 



p»2.«PAI*S2*(2.^f^fMJ2)/XN2 


1)306, 30 1,301 




303 602 = S52+( 2.*wnN) *QQ( I » 


301 If (M2)'>09,3 o9,a^„.„. 

401 C:3I=1,MJ2'''' 'i 

3 PP{I)=0. 

300 3=0* 


K-I-i 




5 

IF ( L 1 "* 1 ) 6 0 » 6 » 6 


60 = 



7 Q( I) =0, 


70 Sl=0. 

Sw * u • . . 

Z2 = l. 

0081 = 1, LI 


Z2»-22^fG( I ) 

^ ‘ ‘ f . 

iu;*^4u >, ; 1 

Sl*SX + {2.#^K)^fQ( I ) 
IIN=M1~I 

do P = 2 . F A I * S I # ( 2 . %J'f ;i! X-Mi 

81 23 = 1 . , . ^ ; 

IF(f.i_l-l)90,9,9 " 

9 0 S S 6 = 0 • 

G0TO12 

y 00101=1, MLl 

10 R(I)=0. 

lOO 23 = 1. _ ,,7 ^ 

- ■ 'rl 

' ^ ' ,,i •; *:■ ,••■• ■■)'■: 








Ak( IT)-CALCRt 


AI (,IT)=*CALCIT 


•■nasi 


( IS. Ml ) .XR( IT 


TO TAfi' £ ON 


dimensional ca 


*'«-MJ2 ) ^fXNl 


IF{J2-1)130, 140,140 


o' LCHECK-Z;,) 144,145. 145, 

144 F=SIN{P) 


GCT0146 


CALE 


(IS) +AR (IT) ^-COS ( P ) -AI’M T 5 *F 


ALC Rb=A 


CALE 


CALC I S = A I ( IS)+AI ( IT)*COS(P)+AR.{lt,l>f 


CALCIS-*^SCALE 


CALCRT = AR ( I $ ) -AR ( I tP ? 



1TC=ISC+N/^ 


R{ ISC)==CALCRS 


FI ( ISC) -CALCIS 


CAlCiT« A I'i.ll ) -1 1 p ) -AR 1 1 T ) 

CALCIT=-CALCIT«SCAU::-' ' 


ARC i3};*<;ALCRS 


lF{J-M 2 )'lla- 5 ?-i'i '52 


152 iSC-ISCi-l 

lTC=ISC+N/2 
FRC ISC)=CALCRS 
F I ( I SC ) =CALC IS 
FRC I TC)=;CALCRT 
_ , ,, FI C I TC)=CALCIT , 

' ' '■ A ^ftYNtSOAtFIxC ISC) ,FI ( ISC) ,FR(nC,) ,FI(ITC) 

IF 

‘ ■‘f . 'R'lf'' 




IS - 


28 3Q( I ) =QQ( I )+l 


IF (CQ( I )-2. ) 304,29,29 


OT028 


iF(L":Ml)30l>3.0l,31. 







'if 

LCHE Ck= LCHECK+l „ t!; 


C'^IFICATION FOR HOLOGRAPHY 


AR( I ) =AFU I )+RR ( I ) 


AR( I ) =AR{ I )**2+AI { I )*#2 


-lODIFICATlON FOR HOLOGRAPHY 




,f QRM A T (A&X»34H60yRieR' 


GC 13600 


XERROR=0 


303d 1=1, N 


AMP( I )=SQRT(D) 


DIPF ( n=C{ I )-AMP( I ) 


R-XERROR+(DlFr( I 


GRT{XERROR/XN) 


5RINT1001 


ORMAl C50X,2 


AI( I ) =0. 





, y. -n/.'J:.' /' 


IN ‘THE >RO0ftAM, 




; ' % , x'V Vj'l' J' 







